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ABSTRACT 

Two  important  results  relating  to  the  uniqueness  of  image  flow  solutions  for  planar 
surfaces  in  motion  are  presented  here.  These  results  relate  to  the  formulation  of  the  im¬ 
age  flow  problem  by  Waxman  and  Ullman  [1],  which  is  based  on  a  kinematic  analysis  of 
the  image  flow  field.  The  first  result  concerns  resolving  the  duality  of  interpretations 
that  are  generally  associated  with  the  instantaneous  image  flow  of  an  evolving  image  se¬ 
quence.  It  is  shown  that  the  interpretation  for  orientation  and  motion  of  planar  surfaces 
is  unique  when  either  two  successive  image  flows  of  one  planar  surface  patch  are  given  or 
one  image  flow  of  two  planar  patches  moving  as  a  rigid  body  is  given.  We  have  proved 
this  by  deriving  explicit  expressions  for  the  evolving  solution  of  an  image  flow  sequence 
with  time.  These  expressions  can  be  used  to  resolve  this  ambiguity  of  interpretation  in 
practical  problems.  The  second  result  is  the  proof  of  uniqueness  for  the  velocity  of  ap¬ 
proach  which  satisfies  the  image  flow  equations  for  planar  surfaces  derived  in-fijCTn  aSP — 
dition,  it  is  shown  that  this  velocity  can  be  computed  as  the  middle  root  of  a  cubic  equa¬ 
tion.  These  two  results  together  suggest  a  new  method  for  solving  the  image  flow  prob¬ 
lem  for  planar  surfaces  in  motion.  s?  ,  ,  ~  / 
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1.  INTRODUCTION 


The  relative  motion  of  a  pin-hole  camera  with  respect  to  a  rigid,  textured 
surface  results  In  a  time-varying  Image  on  Its  projection  screen.  The  Image  flow 
problem  concerns  the  recovery  of  space  motion  and  surface  structure  by  analyz¬ 
ing  Its  evolving  Image  sequence.  Waxman  and  Ullman  [l]  developed  a  formula¬ 
tion  of  the  Image  flow  problem  In  which  the  motion  and  structure  parameters  of 
the  surface  are  given  by  the  solution  of  a  set  of  twelve  non-llnear  algebraic  equa¬ 
tions.  These  equations  were  derived  by  a  kinematic  analysis  of  the  Image  flow 
field.  This  paper  extends  this  earlier  work  of  Waxman  and  Ullman  for  the  special 
case  of  planar  surfaces  In  motion. 

There  are  two  Important  results  presented  In  this  paper  which  relate  to  the 
uniqueness  of  solutions  for  planar  surfaces  In  motion.  Our  first  result  concerns  the 
uniqueness  of  Interpretation  for  orientation  and  motion.  Given  two  successive 
Image  frames  (or  the  Instantaneous  flow)  of  an  Image  sequence.  It  has  been  shown 
by  Waxman  and  Ullman  (l]  that  there  are  generally  two  Interpretations  which 
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are  duals  of  each  other.  In  this  paper  we  show  that  this  duality  can  be  resolved 
to  give  a  unique  Interpretation  when  either  three  successive  Image  frames  (or  two 
successive  flows)  of  one  planar  surface  patch  are  given,  or  two  frames  (or  one  flow 
field)  with  two  planar  patches  moving  as  a  rigid  body  are  given.  A  similar  result 
has  been  obtained  from  a  different  approach  In  the  4-polnt  formulation  of  Tsai 
and  Huang  [3].  Our  approach  to  the  first  case  consists  of  deriving  explicit  expres¬ 
sions  for  the  evolving  solution  over  successive  Image  frames  for  the  two  dual 
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Interpretations,  and  then  comparing  these  expressions  at  a  later  time  when  they 
are  no  longer  dual.  These  expressions  can  be  used  In  practical  problems  to  resolve 
the  ambiguity  In  Interpretation.  In  the  second  case,  we  simply  note  that  the 
spurious  solutions  for  orientations  are  the  same  and  the  rotation  parameters  are 
(In  general)  different  for  all  planar  surfaces  which  are  stationary  with  respect  to 
each  other  In  space.  Our  second  main  result  concerns  the  uniqueness  of  the 
translational  velocity  component  along  the  line  of  sight.  Waxman  and  Ullman's 
solution  of  the  Image  flow  equations  suggested  that  this  velocity  was.  Indeed, 
uniquely  determined  (with  the  duality  described  above  being  the  sole  form  of 
non-uniqueness);  however,  this  was  not  proved.  Here,  we  show  that  given  only  the 
Instantaneous  Image  flow,  the  velocity  component  along  the  line  of  sight  Is 
uniquely  determined.  Further,  It  Is  shown  that  this  velocity  component  can  be 
computed  as  the  middle  root  of  a  cubic  equation.  These  results  are  derived 
directly  from  the  formulation  of  the  Image  flow  problem  In  [l].  Reference  to  a 
similar  result  recently  obtained  by  Longuet-Hlgglns  for  a  different  formulation  of 
the  problem  Is  found  In  Buxton  (4). 

The  organization  of  this  paper  Is  as  follows.  In  the  next  section  we  briefly 
summarize  the  formulation  of  the  Image  flow  problem  of  Waxman  and  Ullman  [i] 
for  planar  surfaces.  In  Section  3,  we  discuss  and  Illustrate  the  dual  nature  of 
planar  surface  solutions.  In  Section  4,  uniqueness  of  Interpretation  for  orientation 
and  motion  of  planar  surfaces  two  Image  flows  Is  proved.  In  Section  5,  the 
uniqueness  of  the  translational  velocity  component  along  the  line  of  sight  for  the 


2.  FORMULATION  OF  THE  IMAGE  FLOW  PROBLEM 

The  relationship  between  a  rigid  object’s  structure  and  motion  and  the  gen¬ 
erated  image  flow  field  has  been  discussed  by  Waxman  and  Ullman(l).  Their  for¬ 
mulation  In  terms  of  “observables”  or  "deformation  parameters”  [2]  Is  central  to 
the  approach  taken  here. 

2.1  Coordinate  System  and  Notation 

We  attribute  the  relative  rigid  body  motion  to  an  observer  represented  by 
the  spatial  coordinate  system  (X,  Y,  Z  )  In  Figure  l.  The  origin  of  this  system  Is 
located  at  the  vertex  of  perspective  projection,  and  the  Z  -axis  Is  directed  along 
the  center  of  the  Instantaneous  field  of  view.  The  Instantaneous  rigid  body 
motion  of  this  coordinate  system  Is  specified  In  terms  of  the  translational  velocity 
V  =  (Vx ,  Vy ,  Vz )  of  Its  origin  and  Its  rotational  velocity 
fl  =  (fijf ,  n  y ,  )•  The  2-D  Image  sequence  Is  created  by  the  perspective  pro¬ 
jection  of  the  object  onto  a  planar  screen  oriented  normal  to  the  Z -axis.  The  ori¬ 
gin  of  the  Image  coordinate  system  (x ,  y )  on  the  screen  Is  located  In  space  at 
(X ,  Y ,  Z )  =  (0,  0,  l);  that  Is,  the  Image  Is  reinverted  and  scaled  to  a  focal 
length  of  unity. 

Due  to  the  observer’s  motion,  a  point  P  In  space  (located  by  position  vector 
R )  moves  with  a  relative  velocity  U  =  -  (V  +  OX  R ).  At  each  Instant,  point 
P  projects  onto  the  screen  as  point  p  with  coordinates 

(x,y)  =  (X/ Z,  Y/Z).  (l) 

The  corresponding  Image  velocities  of  point  p  are  (vx  ,  v )  =  (x ,  y ),  obtained 
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by  differentiating  the  Image  coordinates  with  respect  to  time  and  utilizing  the 
components  of  U  for  the  time  derivatives  of  the  spatial  coordinates  of  P .  The 
result  Is 

[  Vz  Vx  \ 

v z  =  |  x  — - —  j  +  [xy  Qx  -  (l  +  x2)  nr  +  y  nz  }  ,  (2a) 

r  vz  vY  \ 

vy  =  j  y  — - y  J  +  K1  +  y2)  nx  -  *y  nY  -  x  nz  1  •  (2b) 

These  equations  define  an  Instantaneous  image  flow  field,  assigning  a  unique 
2-D  Image  velocity  v  to  each  direction  (x ,  y )  In  the  observer’s  field  of  view.  For 
our  formulation,  we  consider  only  a  single  surface  patch  of  some  object  In  the 

O 

field  of  view. 

2.2  Image  Flow  Equations  for  Planar  Surfaces  in  Motion 

A  small  but  finite  surface  patch  may  be  locally  approximated  by  a  quadric 
surface  described  by  six  parameters:  Its  distance  along  the  line  of  sight  ZQ 
(assumed  to  be  greater  than  zero  throughout),  Its  two  slopes  Tx ,  TY  and  three 
curvatures.  If  this  surface  patch  Is  described  In  our  viewer-centered  spatial  coordi¬ 
nate  system  by  Z  =f(X,Y),  then  using  (1),  Its  local  representation 
Z  =  f(x,y)  can  be  obtained  as  a  second-order  polynomial  In  terms  of  Image 
coordinates  [ij.  Using  this  representation  of  the  surface  and  performing  a 
kinematic  analysis  of  the  Image  velocity  equations  (2)  In  a  small  neighborhood 
around  the  line  of  sight,  we  can  derive  a  set  of  twelve  non-linear,  algebraic  equa¬ 
tions  relating  the  motion  and  structure  of  the  surface  to  a  set  of  "observables”  or 


“deformation  parameters”.  In  these  equations,  the  distance  Z 0  between  the  sur¬ 
face  and  the  camera  along  the  line  of  sight  always  appears  In  ratio  with  the 
translational  velocity  V  and  therefore  Is  not  recoverable  (from  a  monocular  flow 
field).  The  deformation  parameters  are  linear  combinations  of  partial  derivatives 


(up  to  second  order)  of  the  Image  velocity  field  evaluated  at  the  center  of  the 
field  of  view.  They  describe  the  geometrical  distortion  of  the  Image  In  a  small 
neighborhood  around  the  line  of  sight.  For  planar  surfaces  represented  In  the 


form  Z  —  Z0+TxX  +  Ty  Y ,  since  the  curvatures  vanish,  we  get  the  following 
eight  Independent  equations  and  four  consrtalnts: 


—  ~K  -  fiy  > 


—  ~v«  +  » 


-  n  +  v9  tx  , 

-  v,  +  Vy  Tr  , 


•J  (V,  Tx  +  V,  TY) . 


-  -n*  +  f  (V,  Tx  -  V,  Ty) , 

=  -2  (Oy  +  ^  > 


—  nX  ~  VzTy  - 


Og  -  —  Of  ,  012  -  —  07  , 


O  i  o  —  2  Og,  Ou - ~  O  g  . 


where 


(4a,  b) 


(4c, d) 
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for  Zn>  0 


2  \  dxdy  dx2  > 


2  K  dy2  dxdy 


The  quantities  O  x  through  0 12  are  the  Image  deformation  parameters  com¬ 
puted  at  the  origin  of  the  Image  plane.  These  deformation  parameters  can  be 
computed  If  the  Image  velocities  at  four  points  (no  three  of  them  colllnear)  on  the 
Image  plane  are  known  (by  fitting  a  second-order  polynomial  In  Image  coordinates 
to  the  velocity  field;  see  the  next  section.).  Also,  two  methods  of  extracting  these 
deformation  parameters  using  contours  from  successive  Image  frames  have  been 
described  by  Waxman  and  Wohn  [2]. 


The  Image  flow  equations  (3)  form  a  set  of  eight  coupled,  non-llnear  alge¬ 
braic  equations  among  eight  unknowns.  The  method  used  In  [l]  to  solve  these 
eight  equations  Involved  numerically  computing  a  transform  angle  which  aligned 
one  of  the  Image  axes  with  the  direction  of  zero  slope  at  the  point  on  the  surface 
Intersected  by  the  line  of  sight.  In  Section  5  we  describe  a  new  approach  to  solve 
these  eight  equations  which  requires  only  the  solution  of  a  cubic  equation. 


3.  DUALITY  OF  PLANAR  SURFACE  SOLUTIONS 


The  fact  that  the  Image  flow  equations  (3)  are  non-linear  Implies  the  possi¬ 
bility  of  multiple  solutions  corresponding  to  more  than  one  Interpretation  of  the 
scene.  In  fact.  It  was  shown  In  [l]  that  for  any  non-zero  value  of  Vz  satisfying 
equations  (3),  there  are  two  solutions  {VJZ0  ,  O,  Tx,  Ty),  {V/Z 0 ,  O,  Tx,  Ty) 
which  satisfy  the  following  duality  relationship: 


tx  =  -K  /  V.  . 

(7a) 

ty=-V,  /  V ,  , 

(7b) 

UH 

(? 

1 

II 

(7c) 

II 

i 

*■< 

(7d) 

III 

(7e) 

fix  —  nX  ~  Vy  ~  K  Ty  , 

(7f) 

fiy  =  Qy  +  Vx  +  Vz  Tx  , 

(7g) 

fiz  —  0 Z  +  Vs  Ty  ~  Vy  ^x  ■ 

(7h) 

Notice  how  the  slopes  Tx ,  Ty  and  the  components  of  translation  parallel  to  the 
Image  plane  Vx ,  Vy  play  Interchangeable  roles  In  the  two  solutions.  There  are 
two  exceptions  to  this  duality.  The  first  Is  when  Vz  =0  (l.e.  the  velocity  com¬ 
ponent  along  the  line  of  sight  Is  zero)  and  the  second  Is  when  the  translational 
velocity  through  space  Is  parallel  to  the  surface  normal  (l.e.  Vx  /Vz  =  -Tx  and 
Vy/Vz  =-Ty  )  In  which  case  equations  (7)  degenerate  to  Identities.  Many 
numerical  examples  In  [l]  had  Indicated  that  there  was  only  one  value  of  V z 
which  satisfied  equations  (3),  thereby  suggesting  that  there  are  at  most  two  solu- 
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tlons  to  equations  (3).  A  formal  proof  of  this  Is  given  In  Section  4.  We  Illustrate 
the  nature  of  the  dual  solutions  below. 

3.1  Illustration  of  Dual  Solutions 

In  the  recovery  of  surface  structure  and  3-D  motion  from  Image  flow.  It  has 
been  shown  that  It  Is  sufficient  to  describe  an  Image  flow  as  a  locally  second-order 
flow  field  [1,2].  This  has  Implications  with  regard  to  the  surfaces  which  generate 
the  flow  itself.  For  example,  consider  a  planar  surface  patch 

Z  =  Z o  +  Tx  X  +  Ty  Y  for  Z  Q  >  0  .  (8) 

Using  relations  (1),  this  may  be  expressed  In  terms  of  Image  coordinates  as 

Z  —  Z  Q(\-Tx  x  -Ty  y  )~l .  Substitution  of  this  Into  the  image  velocity  equations 
(3)  yields  the  following  expressions  which  are  In  the  form  of  a  second-order  poly¬ 
nomial  In  the  Image  coordinates: 

vx  =  (xVx-Vx  )  (l-Txx-Ty  y )  +  fix  ~  (l+a:2)  fly  +  y  fl^  j  (9a) 
vy  —  {yYx-Vy  )  (i-Ty£ -Ty  y )  +  ^(l+y 2)  fljf  -  ary  fly  -  z  j  (9b) 

For  planar  surfaces,  such  second-order  flows  are  globally  valid  (for  quadric  sur¬ 
faces  the  flow  can  be  locally  approximated  as  second-order).  The  coefficients  of 
this  second-order  flow  then  determine  the  slopes  and  (scaled)  space  motion  of  the 
planar  surface. 

Since  the  above  two  equations  are  globally  valid,  Image  velocities  at  four 
points  (no  three  of  them  colllnear)  or  more  on  the  plane  allow  a  solution  to  be 


obtained.  Of  course  normal  flow  along  contours  (or  edge  fragments)  may  be  used 
as  well  [2].  We  can  graphically  Illustrate  this  duality  of  solutions  by  reexpressing 
the  slopes  In  terms  of  angles: 

Tx  =  tan  9  and  Ty  —  tan  <t>  for  -  —  <  9,  <t>  <  —  ,  (10) 

2  2 

and  then  searching  the  9-<t>  space  for  solutions.  Substitution  of  values  for 
Tx,  Ty  renders  equations  (9)  linear  In  the  unknown  motion  parameters.  There¬ 
fore.  If  Image  velocities  at  more  than  four  points  are  known  (and  the  data  Is 
noisy),  a  linear  least-square  error  minimization  technique  can  be  used  to  solve  for 
the  corresponding  motion  parameters.  A  typical  plot  of  the  least-square  error  In 
the  9-<f>  space  Is  shown  In  Figure  2.  The  actual  error  values  have  been  thres- 
holded  and  Inverted  so  that  peaks  In  the  figure  represent  error  minima.  Loca¬ 
tions  of  peaks  representing  zero  error  (for  noise-free  data)  give  the  solutions  for  9 
and  <t>.  In  the  figure  we  see  that  there  are  two  such  peaks  corresponding  to  two 
dual  solutions.  In  Figure  3  we  Illustrate  how  the  solution  for  structure  (and  con¬ 
sequently  motion)  becomes  unique  as  the  velocity  component  Vg  along  the  line  of 
sight  approaches  zero.  As  V2  —*Q,  one  of  the  two  peaks  corresponding  to  a  dual 
solution  (from  Figure  2)  moves  towards  the  boundary  of  the  9-<t>  space  and  finally 
“disappears”,  resulting  In  a  unique  Interpretation.  This  Is  consistent  with  the 
duality  relations  (6a, b)  In  which  Tx  and  Ty  approach  Infinity.  In  Figure  4  we 
Illustrate  the  behavior  of  the  dual  solutions  as  the  direction  of  the  translational 
velocity  V  approaches  the  direction  normal  to  the  planar  surface.  In  this  case, 
the  peaks  (from  Figure  2)  gradually  move  towards  each  other  and  finally  merge 
Into  one  giving  again  a  unique  Interpretation. 


The  above  method  can.  In  principle,  be  used  to  solve  the  Image  flow  problem 
for  planar  surfaces.  We  need  to  search  the  9-<t>  space  completely  only  once  at  the 
beglnlng  of  the  sequence;  this  search  can  be  done  efficiently  by  first  carrying  out 
the  search  at  coarse  Intervals  and  then  at  finer  Intervals.  For  subsequent  Images, 
we  can  evolve  the  Initial  solution  to  the  next  time  Instant  using  the  expressions 
derived  In  the  next  section,  and  search  a  small  neighborhood  around  this  new 
predicted  solution.  When  the  Image  velocity  data  is  noisy  (e.g.  as  much  as  20% 
perturbations),  we  have  found  It  very  useful  to  first  filter  the  noise  by  fitting 
second-order  polynomials  to  the  flow  field  by  least-square  error  minimization  over 
all  the  points.  The  search  over  9-<j>  space  Is  then  based  on  this  second-order  flow 


3.2  Extraction  of  Deformation  Parameters  from  Points 


Expanding  equations  (9)  and  expressing  them  as  polynomials  In  image  coor¬ 
dinates,  we  find 


vx(x  ,y)  =  a  +bx  +cy  +dx2+d!xy  and 


vy  (x  ,y  )  =  d  +dy  +c'y  +d y2+dxy 


where 


—  -  V,-  fly, 

=  Va  +  V9TX, 

—  nz  +  ?Y* 

=  -  fJ  y-Vz  Tx , 


—  ~Vy  +  nX- 

=  v,+v9ty, 

—  -  Clz  +  Vj,  Tx , 

=  -  ^X~Vi  ty 


(lia) 

(lib) 


(I2a,b) 
(12c, d) 
(12e,f) 
(I2g,h) 


The  eight  coefficients  a,b,c,d9a!,Vfd9d!  can  be  found  If  vx ,  vy  are 
known  at  four  points  (no  three  points  colllnear)  or  more,  by  solving  the  system  of 
linear  equations  (11).  If  these  velocities  are  known  at  more  than  four  points  (for 
noisy  data),  the  coefficients  are  solved  by  a  linear  least-square  error  minimization 
method.  By  Inspection  of  equations  (3)  and  (12),  we  see  that  the  deformation 
parameters  O  x  -  08  can  be  obtained  from  the  eight  coefficients: 

Ox  =  a  ,  02=  o'  ,  03  =  b  ,  04  =  y  ,  (I3a-d) 

°5  =-J-(c+c/),  06=-l(c/-c),  07  —  2d  ,  Os  =  d.  (I3e-h) 
2  2 

From  these  deformation  parameters,  we  can  Immediately  obtain  the  solution  for 
the  3-D  structure  and  motion  parameters  of  the  planar  surface  In  closed  form  by 
the  methods  developed  In  the  following  sections. 


4.  UNIQUENESS  OF  INTERPRETATION 


The  Image  deformation  parameters  0 1  -  08  Introduced  In  Section  2  may  be 
computed  for  a  planar  surface  In  motion  from  a  pair  of  successive  Image  frames 
either  by  tracking  four  or  more  points  or  by  observing  an  evolving  contour  on  the 
surface  [2].  Utilizing  more  than  two  frames  will  allow  a  more  accurate  determina¬ 
tion  of  these  parameters.  These  Image  deformation  parameters  completely 
describe  the  Instantaneous  Image  flow  field.  But  an  Instantaneous  flow  field  In 
Itself  Is  Inherently  ambiguous  as  It  has  two  Interpretations  for  the  orientation  and 
motion  of  the  planar  surface.  This  ambiguity  In  Interpretation  Is  manifest  In  the 
duality  of  solutions  to  equations  (3)  as  discussed  In  Section  3. 

In  this  section,  we  show  that  the  Interpretation  for  structure  and  motion  of 
planar  surfaces  Is  uniquely  determined  In  two  cases.  The  first  case  Is  when  two 
successive  flow  fields  (or  three  successive  Image  frames)  with  one  planar  surface 
are  given,  and  the  second  Is  when  one  flow  field  (or  two  successive  Image  frames) 
with  two  planar  surface  patches  moving  as  a  rigid  body  Is  given.  Notice  that 
from  three  consecutive  Image  frames  we  can  extract  two  sets  of  Image  deforma¬ 
tion  parameters,  one  set  from  each  pair  of  consecutive  frames.  Each  of  these 
deformation  parameter  sets  describes  the  instantaneous  velocity  field  at  the  mean 
of  the  two  time  Instants  (In  the  limit  as  St—* 0)  of  the  Image  frame  pair  from 
which  It  Is  extracted.  We  will  see  below  that,  although  solving  the  Image  flow 
equations  (3)  for  the  first  set  of  deformation  parameters  may  yield  two  solutions, 
only  one  of  them  will  be  consistent  over  time  with  the  second  set  of  deformation 


parameters. 

4.1  Proof  of  Uniqueness  of  Interpretation: 

One  Planar  Surface 

The  scheme  of  our  proof,  when  two  flow  fields  of  a  single  planar  surface  are 
given.  Is  shown  In  Figure  5.  An  Image  flow  field  F  at  time  t  generated  by  a 
planar  surface  In  motion  has  two  solutions,  say 
5  =  {Vx,  Vy,  Vz,  nx,  Hy,  Qz>  Tx,  Ty)  and  Its  dual  SD  which  Is  related  to 
S  by  relations  (7).  For  these  two  solutions,  assuming  that  the  motion  in  space  is 
constant  over  a  short  time,  the  slope  parameters  change  with  time  due  to  rota¬ 
tion  and  the  scaled  translational  velocity  compoments  change  with  time  due  to 
the  combined  effect  of  translation  and  rotation  (we  shall  relax  this  assumption 
shortly).  As  a  result,  after  a  time  6t  later,  S  and  SD  evolve  Into  two  new  solu¬ 
tions,  say  SE  and  SDE ,  respectively.  (Explicit  expressions  for  SE  and  SDE  are 
derived  In  terms  of  S  In  the  Appendix).  The  proof  of  uniqueness  of  Interpreta¬ 
tion  Is  that,  although  S  and  SD  are  duals  of  each  other,  SE  and  SDE  are  not 
(see  Appendix  for  details).  This  Implies  that  If  F'  Is  the  Instantaneous  flow  field 
at  time  t  +6t ,  then  only  either  SE  or  SDE  will  be  one  of  Its  two  solutions.  If 
SE  Is  a  solution  of  the  Image  flow  field  F '  then  the  actual  motion  and  slopes  of 
the  planar  surface  at  time  t  are  given  by  S . 


4.2  Resolving  Duality  for  the  General  Case: 
One  Planar  Surface 


In  general,  for  the  single  planar  surface  case,  the  space  motion  may  not  be 
constant  over  the  time  period  2 St .  We  can  still  proceed  to  resolve  the  duality  as 
follows.  Let  I  v  1 2  and  J3  represent  three  successive  Image  frames  of  an  Image 
sequence  taken  at  very  short  time  Intervals  St  apart.  Then  compute  the  Image 
flow  fields  (possibly  represented  as  Image  deformation  parameters) 
F  and  F '  corresponding  to  the  Image  frame  pairs  (Iv  I2)  and  ( I2 ,  /3),  respec¬ 
tively.  Then  use  the  two  flow  fields  F  and  F'  to  get  dual  solution  pairs 
(S ,  SD)  and  (S',  SD'),  respectively.  Also,  compute  the  time  evolved  solutions 
SE  and  SDE  corresponding  to  S  and  SD ,  respectively,  using  relations  (A12, 
A13).  Now  If  SE  Is  “closer”  to  one  of  S’  or  SD 1  than  SDE  is  to  either  of  them, 
then  5  Is  the  required  solution  of  F  and  the  solution  “closest”  to  SE  among 
S'  and  SD'  Is  the  required  solution  of  F'  .(For  reasons  that  will  become 
apparent  In  the  next  section,  the  determination  of  "closeness  of  solutions”  may 
be  based  entirely  on  Vg .)  Once  we  have  resolved  the  duality  In  this  manner.  In 
subsequent  time  steps  we  need  only  to  repeat  the  process  of  selecting  the  solution 
which  "agrees  most”  with  our  predicted  solution.  Our  predicted  solutions  will  not 
agree  exactly  with  the  computed  solutions  since  the  time  Interval  St  between 
Image  frames  may  not  be  small  enough,  and  the  motion  In  space  may  not  be 
sufficiently  constant  with  time. 


Our  approach  of  seeking  consistent  3-D  Interpretations  over  time  may  be 
viewed  as  a  kind  of  "prediction  and  verification”.  We  evolve  the  solutions  of 
flow  F  forward  In  time  In  order  to  check  for  the  consistent  solution  among  those 
of  F*.  This  Is  an  application  of  the  principle  of  “temporal  coherence"  In  the  3-D 
domain.  Our  approach  may  be  contrasted  with  that  of  Tsai  and  Huang  [3]  who 
select  three  consecutive  Image  frames  Iv  I2  and  J3  ,  and  solve  for  the  parameters 
of  the  backward,  evolution  I2  -*  1 1  along  with  the  forward  evolution  I2  — ►  J3. 
Duality  Is  resolved  by  requiring  consistency  (equality  of  the  slopes)  among  the 
chosen  backward  solution  and  Its  corresponding  forward  solution. 

4.3  Proof  of  Uniqueness  of  Interpretation: 

Two  Planar  Surfaces 

The  proof  of  uniqueness,  when  a  flow  field  with  two  planar  patches  In  the 
scene  Is  given,  Is  quite  straightforward.  We  assume  that  the  two  planar  surfaces 
are  stationary  in  space  with  respect  to  each  other  (l.e.  they  move  as  a  rigid  body) 
and  neither  of  them  passes  through  the  origin  (l.e.  Z0  ^  0  ).  Suppose  that  solv¬ 
ing  the  flow  field  for  the  structure  and  motion  of  the  two  surfaces  yields  the  dual 
solution  pairs  (S(1\  5D(1))  and  ( S(2\  SDW)  where 

5U)  =  (Kju>,  vwt  |r(D,  nx  (1) ,  nr  (1) ,  nz  w ,  Tx  w ,  Ty  (1) )  (i4a) 

sw  =  ( vy\  VW,  vw,  W ,  nr  w ,  nz  m ,  tx  w  ,  Ty  w )  ,  (i4b) 

and  SD  ^  and  SD ^  are  related  to  and  S ^  by  relations  (7).  Then  we  shall 
show  that  the  solutions  corresponding  to  the  actual  Interpretation  for  the  struc- 


ture  and  motion  of  the  two  surfaces,  say  and  S^2\  are  uniquely  determined. 


Suppose  that  the  two  surfaces  are  not  parallel  (l.e.  they  have  different  orien¬ 
tations);  then  for  S W  and  S ^2)  we  have 

(Tx  <*> ,  Ty  «  )  jk  (Tx  <*> ,  ) .  (15a) 

But  for  the  slopes  of  SD^  and  5Z?(2),  from  relations  (5)  and  (7),  we  always  And 
that 

(tx  (1) ,  tY  (1) )  =  (Tx  (2) ,  f^2) )  (15b) 

since 

(-VM/VMt-VM/yM)  =  (-K»/K»(-V»/V«)  ,  (15c) 

which  reflects  the  rigid  body  motion  of  the  planar  surface  pair.  Furthermore, 
S(1)  and  S(2*  have  the  same  rotation  parameters,  whereas.  In  general,  this  Is  not 
true  with  SD(1)  and  SZ>(2).  Therefore,  SDW  and  5Z>(2)  can  be  ruled  out  as  they 
are  Inconsistent. 

If  the  two  planar  surfaces  are  In  fact  parallel  but  distinct  (l.e. 
Z0(1)  jk  ),  then  5(1)  and  5(2)  still  have  the  same  rotation  parameters 

whereas  SD (1^  and  SD *2)  do  not.  Thus,  again  we  can  rule  out  5D(1)  and  5D(2)  as 
unacceptable. 

Therefore,  for  two  (or  more)  planar  surfaces  moving  as  a  rigid  body  (e.g. 
faces  of  a  polyhedron),  a  unique  3-D  Interpretation  may  be  derived  from  an 
Instantaneous  Image  flow  field.  Of  course,  the  flow  field  Itself  must  first  be  seg- 


mented  Into  analytic  regions  corresponding  to  Individual  planar  surfaces  [5]. 


Tsai  and  Huang  [6)  have  arrived  at  a  similar  conclusion  as  a  consequence  of 
their  general  proof  of  uniqueness  for  seven  points  moving  rigidly  through  space. 
Regarding  uniqueness  for  the  specific  case  of  two  planar  surfaces,  our  proof  Is  far 


simpler  and  more  direct. 


5.  PROOF  OF  UNIQUENESS  FOR  Vz 


For  a  planar  surface  In  motion,  from  two  successive  Image  frames,  we  can 
extract  a  set  of  Image  deformation  parameters  which  uniquely  describe  the  Image 
flow  associated  with  the  Image  frame  pair  [2].  Solutions  of  the  Image  flow  equa¬ 
tions  (3)  for  these  deformation  parameters  correspond  to  the  possible  Interpreta¬ 
tions  associated  with  the  Image  frame  pair.  In  this  section,  we  show  that  the 
relative  velocity  of  approach  Vz  associated  with  such  an  Image  frame  pair  is 
unique.  As  a  consequence  of  this  result,  since  there  Is  only  one  pair  of  dual  solu¬ 
tions  associated  with  each  possible  Vz ,  there  are  at  most  two  Interpretations 
associated  with  two  successive  Image  frames. 


5.1  Proof 


We  obtain  an  equation  In  Vz  by  eliminating  all  other  parameters  from  equa¬ 
tions  (3)  as  follows.  From  equations  (3a, b)  we  have 

fly  —  —O  j  —  Vg  and  fix  =  ^2  •  (I®) 

Substituting  for  Ox  and  Qy  from  above  In  equations  (3g)  and  (3h)  and  rearrang¬ 


ing  terms,  we  And 


where 


Vx  =  Vz  T%  -  A  and  Vy  =  Vz  Ty  -  B  , 


A  =  0 1 - 07  and  B  =  02  ~  Oa. 
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Using  the  above  relations,  we  eliminate  Vx  and  Vv  from  (3c-e)  to  get 


>  j-  %  _ _ ^ 


°3  *  vt  +  ( V,  Tx  -  A)  Tx 
04  =  K  +  (V,  Ty-B)  Ty 
2  05  =  (v,  Ty  -B)  Tx  +(V ,  Tx  -  A)  Ty  . 


(18a] 

(19b] 


We  have  here  three  equations  In  three  unknowns.  Equations  (19a,b)  are  qua¬ 
dratic  In  Tx  and  Ty,  respectively,  with  solutions  given  In  terms  of  Vz  as: 


TX+*  TX-  = 


A±y/A*-4Vz(Vg-03) 


(20a) 


ty+>  ty-  = 


B±y/B*-4V,(V,-04) 


(20b) 


Substituting  for  Tx  and  Ty  from  (20)  Into  (19c)  and  simplifying,  yields 


4  OsV,  +AB  =±Vl^2-4n(n-03)l  (B2-4Fz(F,-04)1  .  (21) 


Squaring  expression  (21)  gives 


V,  (  V,3  +  Cx  V*  +  c2  Vz  +  C3  )  -  0  , 


where 


C 1  —  -  (O 3  +  04) , 


c2  = 


C3  = 


-<TA2+T Bl+°i  -o,ot). 

4  4 

—A2  04  +  -B2  03-  -  Os  AB  . 
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(23a) 


(23b) 


(23c) 


In  equation  (22),  we  disregard  the  case  of  Vz  =0  since  substitution  of  this  Into 
equation  (21)  leads  to  an  Identity  not  Involving  Vz  (in  fact,  when  Vz  =  0,  equa¬ 
tions  (3)  have  a  unique  solution  [1]);  therefore  equation  (22)  reduces  to  a  cubic 


equation  In  Vz .  For  convenience,  let  us  express  the  cubic  equation  In  terms  of  a 


I 


new  parameter  a: 


a3  +  C  x  a2  +  C2  a  +  C z  =  0 


We  can  verify  that  the  three  roots  of  equation  (24)  are 

<*o.  «+.«--  V:'j{v*Tx  +  V3TY  +  V,  ±y/l vjl  V*+ W)  (1  +  rj+  W) }  (25a) 

by  expressing  the  coefficients  Cv  C2,  C3  In  terms  of  VK,  Vy ,  Vt,  Tx  and  TY. 
These  three  roots  can  be  expressed  In  the  following  suggestive  way: 


c*o«  a+,  a_  =  Vz ,  -i-  J  v  X  ±  y/(v  •  i/)  (X  ■  X)  J 


(25b) 


where  u  Is  the  scaled  translational  velocity  vector  with  Its  Z  -component 
reversed,  l.e.  u  —  ( Vx  ,  Vy  ,  -  Vz ),  and  X  Is  a  vector  normal  to  the  planar  surface 
given  by  X  =  (  Tx,  Ty*  -l). 


All  three  roots  (25)  are  real,  but  two  of  them,  a+  and  a_,  are  extraneous 
roots  Introduced  by  the  squaring  of  equation  (21).  When  these  roots  are  not 
equal  to  a0,  back  substitution  Into  expressions  (20)  yields  Imaginary  values  for 
the  two  slopes  (l.e.  the  radlcands  are  negative).  When  both  radlcands  are  simul¬ 
taneously  zero  (l.e.  the  direction  of  translational  velocity  Is  parallel  to  the  planar 
surface  normal),  then  relation  (21)  assigns  a  unique  value  to  Vt  given  by 
Vz  =  -  (AB  )/{aO  5)  which  corresponds  to  a0.  In  general,  the  root  a0  Is  the 
unique  solution  for  V, . 


*  •  •  .  •  t  -  .  4  «  *  V  ‘  »*.  •*.  •*. 


5.2  Solving  for  a0  and  the  Orientation  and  Motion  Parameters 


The  fact  that  the  radlcands  In  equations  (20)  are  less  than  or  equal  to  zero 
for  q+  and  a_  and  greater  than  or  equal  to  zero  for  a0  can  be  used  to  show  that 
a_  <  a0  <  a+.  This  relation  Is  useful  In  selecting  a0  from  the  three  roots  of 
equation  (24)  without  actually  substituting  them  Into  the  radlcands  of  equations 
(20).  For  the  root  a0  ,  the  radlcand  In  equation  (20a)  Is  greater  than  or  equal  to 
zero,  l.e.  , 

A  2  -  4  a0  (a0  —  O  3)  >  0 

— *  a02  -  0  3  a0  -  —  A  2  <  0 
4 

-*■  (a0  ~  A_)  (a0  -  A+)  <  0  (where  A  +  ,  A_  —  (1/2)  (03  ±  \J oj  +A  2  )) 

— >  A  _  <  a0  <  A  +.  (26a) 

Similarly,  for  the  radlcand  In  (20b)  we  find 

B.  <  a0  <  B+  (where  B+,B_  =  (1/2)  (04  ±  \/042  4 -B2  ))  .  (26b) 

Hence 

max  (A  B_)  <  a0  <  min  (A  +,  B  +).  (26c) 

Similarly,  using  the  fact  that  the  two  radlcands  are  less  than  or  equal  to  zero  for 
the  roots  a+  and  a_  ,  we  can  derive  the  following  relations: 

a+  >  max  (A+,  P+)  and  a_  <  min  (A  _,  B_)  .  (27) 

Relations  (26)  and  (27)  Imply  that  a0  Is  the  middle  root  of  the  cubic  equation  (24) 
(l.e.  a_  <  a0  S  a+  )• 
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Utilizing  the  closed  form  solution  for  the  roots  of  a  cubic  equation,  we  And 


the  root  a0  of  equation  (24)  as  follows.  Let 

=  £*  _  £1  and  21  2i  2*  +  £± 

P  3  9  aD  q  27  6  2 

Then,  If  p  =  0, 

a0  =  3V-2q  -  b  / 3  ; 


(28a,b) 


(29a) 


If  ?  =  0  , 

a0  =  -b/  3  ;  (29b) 

otherwise. 


a0  =  2  r  Cos  j  j 
where 


b_ 

3 


(29c) 


r  =  slgn(g  )n/|  p  |  and  =  Cos"1 1  —  j  (30a,b) 

Having  solved  for  Vz  =  a0  ,  we  obtain  the  remaining  parameters  as  follows. 
Substitute  this  value  of  Vz  Into  equation  (21)  and  determine  the  sign  of  the  right 
hand  side  which  satisfies  this  equation.  If  this  Is  then  the  dual  solutions  for 
the  slopes  are  ( Tx+  .  Ty+)  and  ( Tx _  ,Tr_);  if  It  Is  the  dual  solutions  for  the 
slopes  are  (Tx+  ,TYJ  and  ( Tx _  ,TY+).  At  this  point  we  can  solve  for  Vx  and 
Vy  using  relations  (17),  and  fly  and  fly  using  relations  (16).  Finally,  Ciz  Is 
found  from  relation  (3f). 


a.  CONCLUSION 


In  this  paper  we  have  extended  the  previous  work  of  Waxman  &  Ullman  [l] 
related  to  the  Image  flow  problem  for  planar  surfaces  In  motion.  We  have 
described  methods  to  disambiguate  dual  Interpretations  associated  with  Image 
flows  generated  by  planar  surfaces.  We  have  also  obtained  closed  form  solutions 
to  the  Image  flow  equations  for  planar  surfaces.  Our  results  suggest  a  new  way  of 
solving  the  Image  flow  problem  for  planar  surfaces,  which  we  summarize  below. 

Assuming  that  the  Image  deformation  parameters  (which  can  be  obtained  by 
either  tracking  points  or  contours  on  the  Image  plane)  at  two  time  Instants  t  and 
t  +dt  are  known,  the  key  steps  of  this  new  method  for  the  single  surface  case  are 
as  follows: 

a)  Initially,  at  time  t  check  If  the  velocity  of  approach  along  the  line  of  sight 
Is  zero  by  noting  if  the  following  condition  Is  satisfied: 

<3i> 

where  A  and  B  are  as  In  equations  (18).  If  the  above  condition  Is  satisfied 
(or  nearly  so  In  the  presence  of  noise),  then  the  Image  flow  equations  (3)  have 
a  unique  solution  which  can  be  easily  computed  directly  from  the  deforma¬ 
tion  parameters  [l].  If  the  condition  Is  not  satisfied,  then  proceed  with  the 
following  steps. 

b)  Compute  Vt  =  aQ  from  the  deformation  parameters  and  then  obtain  the 


remaining  unknowns  as  described  In  Section  5.  At  this  point.  If  there  are  two 
distinct  (dual)  solutions,  proceed  to  the  next  step;  otherwise,  the  unique  solu¬ 
tion  obtained  gives  the  actual  motion  and  orientation  of  the  planar  surface. 

c)  Repeat  steps  (a)  and  (b)  using  the  the  Image  deformation  parameters  at 
time  t+dt ,  and  using  the  method  discussed  in  Section  4  resolve  the  duality 
to  obtain  a  unique  solution. 

If  there  are  two  (or  more)  planar  surfaces  moving  as  a  rigid  body  In  the 
scene,  then  In  step  (b),  resolve  the  ambiguity  of  Interpretation  by  the  method 


described  In  Section  4.3. 


APPENDIX: 


EVOLUTION  OF  PLANAR  SURFACE  SOLUTIONS 

At  time  Instant  t ,  let  a  planar  surface  In  motion  be  described  by 

Z  =  Z 0-\-TxX +  Ty  Y  (Ai) 

In  the  coordinate  system  shown  In  Figure  l,  where  Z0  Is  the  Z  Intercept  of  the 
plane  and  Tx,  Ty  are  the  X  and  Y  slopes  respectively.  Also,  let  V  and  O  be 
the  relative  translational  and  rotational  velocities  In  space  of  an  observer  with 
respect  to  the  planar  surface.  This  relative  motion  In  space  Is  assumed  constant 
over  short  time  Intervals  (over  three  Image  frames).  We  represent  the  relative 
motion  of  the  observer  and  the  structure  of  the  plane  at  time  t  In  matrix  form 
as: 

V  V  V 

vx  r y  rz 

S  =  fix  Cl  y  Cl  2  (A 2) 

Tx  Ty  Z0 

In  this  matrix,  the  first  row  represents  the  translational  velocity  components 
scaled  by  distance  Z0  ,  the  second  row  the  rotational  velocity  components  and 
the  third  row  the  slopes  of  the  plane  and  the  distance  scale  factor.  The  motion 
and  structure  parameters,  represented  by  5  above,  change  with  time  due  to  the 
relative  motion  between  the  plane  and  the  observer’s  reference  frame.  In  the 
remaining  part  of  this  section  we  derive  expressions  for  these  parameters  after  a 
small  time  Interval  dt  has  elapsed. 
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To  compute  the  change  In  the  structure  parameters  during  the  time  Interval 


dt ,  we  first  take  the  time  derivative  on  either  side  of  equation  (Al)  to  get 


dZ_ 

dt 


dZ  o 
dt 


+  T \ 


,  rp  dY  y. 

+  Ty7F  +  Y 


dTy 

dt 


(A3) 


The  time  derivatives  of  (X,  Y ,  Z)  In  the  above  expression  are  given  by  the  three 
components  of  the  vector  -  ( V  +  n XR)  respectively  which  describes  the  rela¬ 
tive  motion  of  a  point  R  —  (X,  Y,  Z)  with  respect  to  the  reference  frame.  Sub¬ 
stituting  these  components  for  the  derivatives,  and  Z0+Tx  X  +  Ty  Y  for  Z  and 
rearranging  terms,  we  get 


|  ~  ZQ  ( (fly  +  VX/Z0)TX  -  {Px-Vy / Z 0)Ty  -  Vz /Z0  )j  + 

|  -  (  Tx(fly  TX  -  fly  Ty )  +  (Oy  +>  S\Z  Ty)  )  j  X  + 

{  -  (  TyiCly  TX  -ClXTy)~  (ClX  +  UZ  TX)  )J  Y  =  0. 


In  the  above  expression,  since  X ,  Y  are  Independent  parameters  of  points  on  the 
plane  In  motion,  we  can  equate  each  of  the  three  terms  above  to  zero  separately 
to  get  the  exact  differentials  for  the  slopes  and  distance  as 


II 

o 

$ 

o  £  (fiy  +  V,  )TX  -  (nx-Vy)TY  -  Vz  j 

dt  , 

(A4a) 

dTx  = 

Ty  ~  &X  Ty)  +  (Oy  +  n  z  Ty, 

)]  dt  , 

(A4b) 

dTy  = 

Ty(QyTX  -  Qx  Ty)  -  (H*  +  ClZ  TX  ) 

j  dt  . 

(A4c) 

Using  the  above  relations,  we  can  compute  the  new  structure  parameters  at  time 

t  +dt  as 


TJc  =  Tx  +  dTx  ,  Tfr  =  Ty  +  dTy  and  Zq  —  Zq  +  dZ0  .  (A5a-c) 
The  new  translational  velocity  V'  at  time  t+dt  Is  found  (In  the  absence  of 


accelerations)  from 

V>  =  V  +  Vxn  dt  .  (A6) 

Dividing  V'  by  Zq  and  simplifying,  we  get  the  new  scaled  velocity  components 
Vi  =  vx  +  (V9  Clz  -  VznY  -  vx  s)  dt  +  O  (dt2) ,  (A7a) 

Vy  =  Vy  +  ( VZ  &X  -  VX  &Z  ~  Vy  s)  dt  +  O  (dt2) ,  (A7b) 

Vi  =  vt  +  (vx nY  -  v3 nx  -  vt  s)dt  +  o (dt2) ,  (A7c) 

where 

$  =  (flY  +  V,  )TX  -  (Hy  -  Vy  )Ty  -  Vz  .  (A8) 

The  rotational  velocity  remains  unchanged  since 

n'  =  n  +  nxn  dt  =  n  .  (aio) 


Summarizing  the  above  results,  we  have  that,  if  the  structure  of  a  planar 
surface  and  the  Instantaneous  motion  of  an  observer  relative  to  the  planar  surface 
at  time  t  Is  S  given  by  (A2),  and  at  time  t  +dt  Is  SE  given  by 


SE  = 


Vi  Vi  Vi 


n  n  rj 


then  S  and  SE  are  related  as  below. 


V 

VM 

vi 

= 

Vy 

+ 

Vi 

vt 

— s  n  2  "*n  Y 

-ci2  -s  nx 

O  y  -tljf  —S 


dt  , 


28 


(All) 


(Al2a) 


(A12b) 


(  rfy-,  f2 y,  Qz  )  =  (  Hy*  Qy,  nz  )  and 

Tfr  Tx  nYTx-nxTY  nz  nY 

Tj-  =  Ty  +  -nz  ^yTx-^xTy  ~&x 

Z'o  z0  Z0(ClY  +  Vx)  -Z0(Clx~Vy)  -ZQVz 


Tx 

Ty  dt.  (A12c) 


Now  we  refer  the  reader  to  Figure  5  for  the  scheme  of  our  proof  of  uniqueness. 
From  our  discussion  In  Section  3,  the  Instantaneous  Image  velocity  field 
corresponding  to  S  has  another  solution  SD  which  Is  the  dual  of  5  given  by 


\-TxV, 


-tyv. 


sd  =  n x-v,-vzty  ciy+vx+v,tx  n z  +  vxty-vstx 


(A13) 


-Vt/Vx 


~VJVX 


We  have  expressed  SD  and  SE  In  terms  of  S .  Let  SDE  denote  the  new  solu¬ 
tion  obtained  by  evolving  SD  In  time  to  t+dt.  SDE  can  be  obtained  by  substi¬ 
tuting  the  elements  of  SD  for  the  corresponding  elements  of  S  In  expression 
(A12)  for  SE .  The  proof  of  uniqueness  of  Interpretation  Is  that,  at  time  t , 
although  S  and  SD  are  the  two  dual  Interpretations  of  a  planar  surface  In 
motion,  at  time  t+dt ,  their  time  evolved  solutions,  SE  and  SDE  respectively, 
are  not  duals.  A  simple  way  to  verify  that  SE  and  SDE  are  not  duals  Is  to  note 
that  their  Z -components  of  translational  velocity  are  not  the  same  (for  dual  solu¬ 
tions,  this  should  be  the  same.).  It  can  be  shown  easily  that  the  values  of  this 
velocity  component  for  SE  and  SDE  are  related  by 


v,  (SOS,  -  K  («>  +  (  V,2  -  K‘  TJ  +  VJ  -  V ,*  TJ  )  dt 


(AM) 


The  two  velocity  components  are  the  same  (l.e.  the  term  with  dt  above  vanishes) 


when  the  direction  of  translational  velocity  Is  parallel  to  the  plane's  normal.  In 
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Figure  2.  Illustration  of  least  square  error  In  the  0  -  <j>  space.  The  two  peaks  In 
the  figure  represent  error  minima  and  their  locations  give  the  two  dual  solutions. 
Here,  the  scaled  velocity  component  along  the  line  of  sight  Is  Vz  =  0.4  and  the 
angle  between  V  and  the  surface  normal  Is  0  =  93° . 


(a)  (b) 

Figure  3.  Behavior  of  dual  solutions  as  the  velocity  component  along  the  line  of 
sight  approaches  zero.  Beglnlng  from  Figure  2,  we  see  that,  as  VT  — *  0,  one  of 
the  peaks  moves  towards  the  edge  of  the  0  -  <f>  space  and  eventually  "disappears” 
resulting  In  a  unique  solution,  (a)  Vz  =  0.1  .  (b)  V z  =  0.01  . 


(a)  (b) 

Figure  -4.  Behavior  of  dual  solutions  as  the  direction  of  V  approaches  that  of  the 
surface  normal.  Beglnlng  from  Figure  2,  we  see  that,  as  (3  — ►  0,  one  of  the  peaks 
moves  towards  the  other  and  eventually  merges  with  It  resulting  In  a  unique  solu¬ 
tion.  (a)  /?  =  67°  .  (b)  /?  =  0°  . 
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Figure  5.  The  uniqueness  of  the  3-D  Interpretation  for  planar  surfaces  In  motion 
stems  from  coherence  over  time. 
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